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Abstract. Calculations in a (3+l)-dimcnsional model indicate that Pauli- Villars reg- 
ularization can be combined with discrete light-cone quantization (DLCQ) to solve at 
least some field theories nonperturbatively. Discrete momentum states of Pauli- Villars 
particles are included in the Fock basis to automatically generate needed counterterms; 
the resultant increase in basis size is found acceptable. The Lanczos algorithm is used 
to extract the lowest massive eigenstate and eigenvalue of the light-cone Hamiltonian, 
^C) • with basis sizes ranging up to 10.5 million. Each Fock-sector wave function is com- 

puted in this way, and from these one can obtain values for various quantities, such as 
average multiplicities and average momenta of constituents, structure functions, and a 
2^ ' form-factor slope. 

9*: INTRODUCTION 

Ph. 

Field-theoretic calculations of bound-state properties, such as those one would 
like to do for quantum chromodynamics (QCD), require regularization of infinities 
and renormalization of parameters. As a way of providing a systematic regulariza- 
tion of ultraviolet infinities we study [1,2] the use of Pauli- Villars (PV) regulariza- 
tion [3] in the context of discrete light-cone quantization (DLCQ) [4,5]. Renormal- 
ization is accomplished by adjusting bare parameters to fit selected state properties 
with "data." The problem to be solved is then a bound-state eigenvalue problem, 
which includes PV constituents, combined with renormalization conditions. The 
couplings of the PV constituents are chosen to produce desired cancellations in 
perturbation theory. 

We have tested these ideas for two related (3+l)-dimensional Hamiltonians [1,2]. 
The first [1] was constructed to have an analytic solution, in analogy with the 
equal-time model of Greenberg and Schweber [6]. The second Hamiltonian [2] is a 
generalization of the first which assigns proper light-cone energies to all particles, 
but does not have an analytic solution. Both Hamiltonians are distantly related to 
Yukawa theory, in that a fermion field acts as a source and sink for bosons. The 
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second Hamiltonian and some of the results obtained will be described here. Work 
on direct application to Yukawa theory is in progress. 

The choice of light-cone coordinates (x ± = t ± z, xj_ = (x, y)) [7,5] is driven 
by important advantages, which include kinematical boosts, a simple vacuum, and 
well-defined Fock-state expansions with no disconnected pieces. The latter two 
derive from the positivity of the longitudinal light-cone momentum p + = E + p z . 
When momenta are discretized [4] this positivity brings the additional advantage 
of a finite limit on the number of constituents. 



MODEL EIGENVALUE PROBLEM 



The bound-state eigenvalue problem is H LC ®a = M 2 Q a , where H LC = P + P~ is 
known as the light-cone Hamiltonian, P + is the longitudinal momentum operator, 
and P is the generator for evolution in light-cone time. We work in the frame 
with no net transverse momentum and in a basis diagonal in P + . 

The Hamiltonian that we consider is 
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The creation operators al, and associated with fermion, boson, and 

PV boson fields, respectively. The corresponding masses are M, //, and fi\. Each 
operator depends on a light-cone three-momentum such as p = (p + ,p x ,P y )- The 
nonzero commutation relations are 



= 16n 3 q + 5(q - q') , 
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= 167T 3 q + 5(q - q') . 



The structure of the Hamiltonian provides for emission and absorption of bosons 
by the fermion, but no change in fermion number. We explore only the one-fermion 
sector. The particular form of the interaction causes the fermion mass counterterm 
to have an unusual momentum dependence. The coefficient of this counterterm is 
finite because of cancellations arranged by assigning an imaginary coupling to the 
PV boson. 

The state vector describes a dressed fermion with spin a. Its Fock-state 
expansion is given by 
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To satisfy the eigenvalue condition, the Fock-sector wave functions (n ' ril ) must 
solve the following coupled system: 
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The bare parameters M' and 5 are determined by fitting (:0 2 (O):) = $£. :0 2 (O): $ CT 
and M 2 to chosen values. The quantity (:0 2 (O):) was selected for ease of computa- 
tion; it can be computed from a form similar to the normalization sum (4): 
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The renormalization conditions that determine M' Q and 5 are then solved simul- 
taneously with the eigenvalue problem (5). In practice this is done by rearranging 
(5) into an eigenvalue problem for 1/g and simultaneously solving for M' Q in a sin- 
gle nonlinear equation where (:0 2 (O):) is equal to a fixed value. The simultaneous 
solution is done by iterative means. 

Once the wave functions 0( n ' ni ) have been obtained, they can be used to compute 
various quantities, such as the boson structure function 
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and the average boson multiplicity and momentum 

N = / fB(y)dy, (y) = f yf B {y)dy. 
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We also compute the slope of the fermion "charge" form factor -F'(O), from an 
expression derived in Ref. [1]. 



NUMERICAL METHODS AND RESULTS 

The coupled equations (5) are converted to a finite matrix eigenvalue problem 
by applying the DLCQ procedure [4]. Integrals are approximated by sums over 
discrete momentum values (nn/L, n x 7i/L±, n y ir/L±), and the transverse range is 
limited by a cutoff A 2 such that (mf + p±i) < A 2 p~l / P + for each particle, with 
rrii being its mass. The length scales L and L± are associated with a light-cone 
coordinate box. Bosons are assigned periodic boundary conditions and the fermion 
is assigned an antiperiodic boundary condition in the longitudinal direction. The 
momentum integer n is then even for bosons and odd for the fermion. 

The total longitudinal momentum P + of the dressed fermion defines an odd 
integer K = P + L/tt, called the harmonic resolution [4]. A longitudinal momentum 
fraction x = p + /P + then reduces to a rational number n/K. The positivity of 
longitudinal momentum implies that < n < K and that the maximum number of 
constituents is of order K/2. The integers n x and n y range between — N± and N±, 
with N± set to reach the limit imposed by the transverse cutoff for the one-boson 
physical state. 

Thus the discretization is determined by three parameters: K, Nj_, and A 2 . 
The transverse scale L± is computed from these as L± = itN±\J2/ (A 2 — M 2 — ji 2 ). 
The longitudinal scale does not appear [4], but the limit L — > oo is equivalent to 
K — * oo. We therefore study the limit where K, N±, and A 2 all become large. This 
recovers the continuum form of the theory, which is regulated by the PV mass fj,i, 
not by A. We must then also study the large Hi limit. 

Typical discretizations, such as K = 15, N± = 6, and A 2 = 50/i 2 , with /i 2 = 10fi 2 
produce matrices with ranks on the order of 5 million. The largest calculations 
carried out were of rank 10.5 million, which required approximately 2 hours of cpu 
time on a single 4-processor node of an IBM SP. The diagonalization method used 
is the Lanczos algorithm for complex symmetric matrices [8]. 

Although the automatic truncation of particle number imposed by DLCQ can be 
sufficient, further truncation can be made when the coupling is weak. Such trun- 
cation, typically to 4 bosons, was used to permit increased resolution within fixed 
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memory limits. The validity of such an approximation was checked by computing 
the contribution of individual Fock sectors to the total norm. 

Most quantities are remarkably insensitive to numerical resolution. This can be 
seen in Fig. 1 where we display the boson structure function for various mass 
ratios Mj fi. Different values of K and N± generally do not yield significantly 
different results. Notice that a smaller mass ratio is associated with a state where 
the boson constituents carry more of the momentum. 

Values for a set of bound-state observables are given in Table 1. The quantity 
1 0o 1 2 represents the probability for the bare fermion state. Each entry has been 
extrapolated from the numerical results by fits to the form a + (3/K 2 + 7/-/V 2 . To 
obtain the behavior of this form the use of weighting factors [1] in DLCQ sums is 
important. The numerical resolutions used ranged from 9 to 19 for K and 5 to as 
much as 10 for N±; the larger values of Nj_ were available only for smaller K. Most 
observables converge quickly with respect to the PV regulator mass Only M' Q 
is strongly dependent on The form factor slope is sensitive to the transverse 
resolution and range. 

SUMMARY 

This work shows that PV regularization is feasible for DLCQ calculations. The 
matrix size does increase but not beyond the capacity of present-day machines for 
the models considered. More complicated theories will require multiple computing 

0.30 n 1 




FIGURE 1. The boson structure function fs at various numerical resolutions and mass values, 
with (:(f) 2 (0):) = 1, A 2 = 50/j 2 , and [i\ = 10^ 2 . The solid lines are parameterized fits of the form 
Ay a (l - y) b e- c y. 
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TABLE 1. Extrapolated bare parameters and observables, with (:0 2 (O):) = 1. 



(M//z) 2 

W^) 2 
(a/m) 2 


1 

5 

12.5 


1 

5 

25 


1 

5 

50 


1 

10 

25 


1 

10 
50 


1 

20 
50 


1 

20 
100 


0.1 
10 
50 


5 

10 
100 


10 
10 
100 


M'/u 2 


21.4 
1.26 


17.7 
1.10 


16.3 
1.10 


17.8 
1.48 


16.0 
1.4 


16.0 
1.8 


15.5 
1.9 


15.1 
1.39 


18.1 
1.66 


19.0 
1.60 


\<f>0? 

-l00u 2 F'(0) 


0.82 
1.04 
0.18 


0.83 
0.78 
0.15 


0.84 
0.66 
0.14 


0.85 
0.72 
0.15 


0.86 
0.59 
0.14 


0.87 
0.59 
0.13 


0.87 
0.51 
0.13 


0.83 
2.0 
0.16 


0.89 
0.14 
0.10 


0.90 
0.07 
0.09 



(y) 0.077 0.062 0.057 0.062 0.056 0.056 0.053 0.073 0.032 0.024 

nodes and message-passing technology. Use of multiple nodes is facilitated by a 
natural block structure that arises in the matrix due to limited coupling between 
Fock sectors. 

Work on Yukawa theory in a single-fermion truncation is now in progress. The 
complications include additional PV boson flavors and nontrivial spin dependence. 
Quantum electrodynamics is perhaps the next logical step. QCD could also be 
considered in a broken supersymmetric form that contains heavy particles analogous 
to the Abelian PV particles introduced here. 
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